Low-dimensional weakly interacting Bose gases: non-universal equations of state 



0^ 
O 
O 
(N 



O 



C/2 

cr: 



I 

o 



> 
m 
in 

o 

ON 

o 



■ 

X 



G.E. Astrakharchik^, J. Boronat^, I.L. Kurbakov^, Yu.E. Lozovik^, F.Mazzanti^ 

^ Departament de Fisica i Enginyeria Nuclear, Campus Nord B4--B5, 
Universitat Politecnica de Catalunya, E-08034 Barcelona, Spain and 
^ Institute of Spectroscopy, 14-2190 Troitsk, Moscow region, Russia 
(Dated: July 30, 2009) 

The zero-temperature equation of state is analyzed in low-dimensional bosonic systems. In the 
dilute regime the equation of state is universal in terms of the gas parameter, i.e. it is the same 
for different potentials with the same value of the s-wave scattering length. Series expansions of 
the universal equation of state are reported for one- and two- dimensional systems. We propose 
to use the concept of energy-dependent s-wave scattering length for obtaining estimations of non- 
universal terms in the energy expansion. We test this approach by making a comparison to exactly 
solvable one-dimensional problems and find that the generated terms have the correct structure. 
The applicability to two-dimensional systems is analyzed by comparing with results of Monte Carlo 
simulations. The prediction for the non-universal behavior is qualitatively correct and the densities, 
at which the deviations from the universal equation of state become visible, are estimated properly. 
Finally, the possibility of observing the non-universal terms in experiments with trapped gases is 
also discussed. 



I. INTRODUCTION 

Understanding the properties of rarefied quantum sys- 
tems is a fundamental question that has been addressed 
in a large number of works. This problem was exten- 
sively studied in the 50s-60s when significant develop- 
ment of mathematical formalism (perturbative methods, 
Feynman diagrams, diagonalization techniques, etc., see, 
for example, [EH]) permitted to obtain important re- 
sults and brought a lot of interest to dilute quantum sys- 
tems. Some important results were as well obtained in 
low-dimensional systems |^, which at that moment 
were rather mathematical toys with reduced applicabil- 
ity in the real world. The situation changed radically 
with the realization of Bose-Einstein condensation in di- 
lute gases [5i,J3j. Having an excellent experimental control 
over the geometry of the cloud it was possible to create es- 
sentially pure quantum gases in the dilute regime and to 
probe the system properties. The experimental advances 
in the field with the realization of very anisotropic traps 
stimulated further the interest in dilute low dimensional 
gases (see, for example, [7l.[8l.l9l.[Tol.[TTI.[l^). 

In the ultradilute limit the interparticle potential can 
be described by one parameter, namely the s-wave scat- 
tering length a, and the ground state properties of a gas 
are governed by the gas parameter na^ , where n is the 
particle density and D stays for the dimensionality. As 
the density is increased details of the interaction poten- 
tial become important. Such a non- universal regime has 
been thoroughly studied in three-dimension al g eometries, 
where the universal terms are known [l^, Il4l |. Low en- 
ergy corrections coming from the specific interaction po- 
tential can be described by the effective range tq. Cor- 
rections to the ground-state energy, excitation spectrum 
and condensate fraction can be obtained ( see, for exam- 
ple, [H, [l^ and more recent works [13, [i3|)- It was 
shown that for two body problem the inclusion of an 
energy-dependent pseudopotential improves significantly 



upon the use of an energy-independent pseudopotential 
Also, the concept of momentum dependent scatter- 
ing length is very useful for estimation of the interaction 
for a Rydberg atom where it allows to take into account 
the effect of the Coulomb potential of nucleus, see for ex- 
ample Ref f20l| . Unfortunately, much less is known in low- 
dimensional systems. Indeed, only recently the universal 
terms of the 2D equation of state have been correctly 
derived [111, m, nil and checked numerically [l^]. In the 
present study we address the problem of non-universal 
corrections in low-dimensional systems. 



The rest of the article is organized as follows. In 
Sec. |TT] we discuss the origins of the universal behavior 
and study the two-body scattering problem and propose 
a simple way to obtain non-universal corrections. In 
Sec. mil the equation of state of some exactly solvable 
one-dimensional models are analyzed. Some properties 
of two-dimensional systems are addressed in Sec. IIVI We 
start with an overview of the literature in Sec. IIV Al In 
Sec. lIVBl we discuss the expansion of the universal equa- 
tion of state and provide some physical insight on the ori- 
gins of the beyond mean- field (BMF) terms. The knowl- 
edge of the expansion of the universal equation of state 
permits us to investigate the non-universal equation of 
state as it comes from the method proposed in Sec. HI] 
and confront that with numerical results. Section IIV CI 
is devoted to the study of non-universal effects in the s- 
wave scattering problem and in the many-body equation 
of state. In Sec. |V] we discuss the possibility of exper- 
imental observation of non-universal effects in trapped 
cold gases. The feasibility of reaching an ultradilute two- 
dimensional regime is also discussed. Finally, the main 
conclusions are drawn in Sec. IVII 
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II. UNIVERSAL AND NON-UNIVERSAL 
TERMS 

In dilute systems the probability of three-body colli- 
sions is highly reduced leaving the two-body scattering 
the most important physical process. In this process two 
particles scatter each other with a relative momentum 
k. The two-body scattering problem is described by the 
Schrodinger equation 
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where /i is the reduced mass. If the interaction potential 
Vintif) is short-ranged, its exact shape is not important 
at low density and the relevant quantity of the scattering 
solution tlj{r) is the phase 6{k) at distances larger than 
the range of the potential. For small scattering energies 
the phase can be expanded in terms of the momentum k. 
In a 3D system this leads to 



fcctg(5(fc) = -— + ifc^ro 



(2) 



where ao is the s-wave scattering length and tq is the 
effective range. If the scattering momentum is very 
small the only relevant parameter is the s-wave scatter- 
ing length ao and all potentials having the same value of 
ao will behave similarly. This limit is known as universal 
regime. The relevant length scales are then ao and the 
interparticle distance. It is expected that the many-body 
ground-state energy can be expressed in terms of the gas 
parameter naQ , where D denotes the dimensionality of 
the problem. 

For example, the low density energy per particle of a 
homogeneous weakly-interacting Bose gas in 3D at zero 
temperature is given by 
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with the leading term linear in the density being the 
mean-field Gross-Pitaevskii contribution 25] and quan- 
tum fluctuations contributing to the subleading n'^/^ Lee- 
Huang- Yang correction[l^[l3. The next term scales like 
n?, but it is no longer universal [Tsl [iGj and depends on 
the explicit choice of the interaction potential. 

It is possible to recast the definition ^ of the scatter- 
ing length ao in a different form, namely, as the position 
of the node of the analytic continuation of the scattering 
solution from distances much larger than the range of the 
potential in the zero-energy scattering limit. Indeed, in 
3D, in the limit of very low-scattering energy the phase 
reads 5{k) = —kao and the scattering solution becomes 
sin(A:r + 6)/r k{r — cLo)/r, which has a node at r = ao. 
The advantage of the alternative definition is that it is 
well suited also to low-dimensional problems. 

We generalize the last definition to finite values of the 
scattering energy. 



Definition 1. The generalized scattering length a{k) is 
the position of the node of the analytical continuation of 
the large distance r — > cx3 two-body scattering solution 
'ip{r) at the scattering energy h^k^ /2^. If there are sev- 
eral nodes, the position of the closest node to r = is 
considered. 

In this way the s-wave scattering length a{k) depends 
on the scattering momentum and fulfills the condition 
limfc^o = o-Q. An example how the s-wave scatter- 
ing length changes with the type of the potential is shown 
in Fig. [1] for several characteristic interactions in one di- 
mension. The figure shows the asymptotic continuation 
of zero energy scattering solution. We will show in the 
next sections that the inclusion of the finite-momentum 
corrections improves the description of the energy and 
allows us to estimate the term of the expansion where 
the non- universal behavior appears. 

At this point it is important to understand the re- 
lation between the effective range and the s-wave scat- 
tering length in the description of non-universal effects. 
The effective-range theory is well established in three- 
dimensional systems (see, for example, textbook [T]). 
The effective range is then defined from the expansion of 
the phase shift in terms of the scattering momentum, see 
Eq. ([2|) . The constant term defines the s-wave scattering 
length ao, and the effective range vq corresponds to tak- 
ing into account dependence on k^. Instead, the energy 
dependent s-wave scattering length a{k) includes in ad- 
dition all higher order momenta, i.e. fc^, fc*, k^,... More 
importantly the concept of a{k) can be applied to low- 
dimensional systems, where the non-universal terms in 
the equation of state are not generally known. In our ap- 
proach it is enough to know the dependence on ao of the 
universal equation of state and the non-universal terms 
will be automatically generated. 



III. ONE-DIMENSIONAL SYSTEMS 

One peculiarity of the one-dimensional world is that 
several many-body models can be solved exactly (with 
short- [1, |3| and long- [2^ range interactions), in the 
sense that the exact ground state can be written either 
exphcitly[3L[2^ or can be easily obtained as the solution 
of a system of integral equations This allows us to 
test the proposed approach of using an energy dependent 
scattering length by comparing to the exactly known re- 
sults. 

The ground-state energy of a Bose gas with a repul- 
sive (5-pseudopotential interaction (Lieb-Liniger model) 
can be obtained by solving Bethe ansatz equations. The 
expansion of the energy in the mean-field regime Q has 
a structure similar to that of the three-dimensional case 
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where giD = —2h'^/{maiD) > is the one-dimensional 
coupling constant. Indeed, the leading term in Eq. ^ 
is the same as it would come out from the mean-field 
Gross-Pitaevskii equation, while the subleading term is 
the same as obtained from Bogoliubov theory. In pass- 
ing by we note that such a coincidence is not obvious a 
priori, as both the Gross-Pitaevskii and Bogoliubov the- 
ories assume that all or a large fraction of particles are in 
the condensate. Instead, strictly speaking, Bose-Einstein 
condensation in homogeneous one-dimensional system is 
absent [27). 

The reason why the theories based on the presence of 
a Bose condensate produce correct results for energetic 
properties can be understood by following the similar 
arguments used in the renormalization group approach 
(see, e.g. Ref. [2lj). The main contribution to the energy 
comes from short distances. At short distances the phase 
coherence may be present even in the absence of the Bose- 
Einstein condensation. Therefore, on this length scale 
it is possible to apply the perturbative theories that are 
based on the assumptions of a macroscopic occupation of 
the condensate. Coherence at finite distances larger than 
the interparticle distance is sufficient for MP and Bogoli- 
ubov theories to yield correct result for the ground-state 
energy. In particular, such theories successfully describe 
one-dimensional systems at zero temperature (such as 
Lieb-Liniger gas in the regime of weak correlations) and 
two-dimensional dilute Bose gas at finite temperature, 
none of which has true Bose-Einstein condensation. A 
mathematical way to resolve the paradox and to prove 
the validity of the Bogoliubov result in one-dimensional 
systems is to use space discretization and to introduce 
the concept of a quasi-condensate|29j. 

Contrary to three- and two- dimensional systems, here 
the mean- field regime means high densities nioaio 3> 1. 
This precludes us from using the concept of energy- 
dependent scattering length in the MF regime, as the 
energy of an incident particle would be huge, see Eq. ([5]), 
and so would be the deviations of a{k) from oq. Thus, 
the mean-field regime is no longer universal (contrary 
to what happens in 3D and 2D systems), as the energy 
and correlation functions are very different for the 6- 
pseudopotential [l,[30,[3l|, the Calogero-Sutherland 
potential [1^, H^l and the dipolar l/\z\^^ interaction [s^ j. 
Instead, in the regime of strong quantum correlations, 
?T.|ai£)| <C 1, the energy and correlation functions of all 
those models (essentially, for any repulsive interaction 
potential) approach the same universal limit referred as 
Tonks- Girardeau [sj regime (see also Fig. [I]). 

This is a peculiarity of the one-dimensional world that 
the dilute regime, ?T.i_D|ai£)| <^ 1 corresponds not to a 
mean- field limit, but rather to a regime where quan- 
tum fluctuations are dominant. The energy in this limit 
is given by the energy of an ideal Fermi gas E/N = 
TT^h^n^j2,/{6m) and the wave function of strongly inter- 
acting bosons can be mapped onto a wave function of 
non-interacting fermionsfsl. |34|. [35| . For instance, the en- 
ergy of a gas of hard rods of size am > is obtained from 



the energy of an ideal Fermi gas by taking into account 
the excluded volumejl]: uid — > ^/{L — Nam)- We ex- 
pand this expression in terms of the one-dimensional gas 
parameter h — ni^aiD at small densities n <C 1 and get 
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It is worth mentioning that the beyond-mean-field 
terms in three-dimensional systems were first obtained 
for a hard-sphere gas by Lee, Huang, Yang^^, and 
afterwards were shown to be universal [1, [13, 3|- Start- 
ing from the expansion for a one-dimensional analogue 
for hard-spheres, Eq. ((S]), we will calculate the first non- 
universal corrections for a different potential. We chose 
a (5-pseudopotential, as its exact groundstate energy is 
known and thus we can test our approach. The solu- 
tion of the scattering problem ((!]) with Vintir) = giD^ir) 
can be readily written 'ip{r) oc sin(fc|r| — arctg kaio)- The 
energy-dependent s-wave scattering length can be explic- 
itly expressed as a function of the momentum and the 
leading correction to oq is quadratic in momentum 
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The substitution of ^ into (O for a characteristic value 
of the energy h'^k'^/m cx n'^h'^nl^/{6m) allows us to es- 
timate the first correction due to non-universality: 
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(7) 



A possible concern about the validity of the obtained 
result is that expansion ([5]) is done for tiidOid > 0, 
while expansion ([7]) is used to describe a region where 
niDO-iD < 0, with a different sign of the s-wave scatter- 
ing length. We argue that the universal equation of state 
is smooth as a function of the one-dimensional gas param- 
eter niuaiE). This is supported by the apparent similari- 
ties between the hard-rod gas and the gas- like state of the 
attractive (5-pseudopotential ( "super- Tonks-Girardeau" 
system) [s^ . We also note that the Bethe ansatz solution 
for two-component attractive and repulsive fcrmions is 
continuous (compare results of Refs. [37, 38, 39]). 

The result can be compared to the exact predictions 
for the Lieb-Liniger model based on the Bethe ansatz 
technique. The exact result can be obtained by solving 
the integral equations recursively (details are given in 
Appendix \K\ and reads 
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By comparing the exact results for the hard- rod gas, 
Eq. (O, exact results for the (5-pseudopotential gas, 
Eq. and the approximate result ([7]) obtained by the 
proposed method we conclude that: 

• The order of the expansion in which the 6- 
pseudopotential and hard rod energies differ is pre- 
dicted correctly 
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Figure 1: Solid lines: typical two-body scattering solutions 
^(r) at zero energy for Lieb-Liniger (upper curve), Tonks- 
Girardeau (middle curve) and Hard-Rod (lower curve) Hamil- 
tonians. Dashed line: analytic continuation of the scattering 
solution for the Lieb-Liniger model. Arrows point the posi- 
tions of the nodes. 



• The expansion ([8]) contains same rational terms as 
the expansion ([5]), while in addition it has irrational 
terms (here multiples of tt^). The use of energy- 
dependent scattering length permits to guess cor- 
rectly the structure of the potential-dependent cor- 
rection. 



We find that the first three terms of the expansion 
are the same for considered potentials. The physical 
meaning of such terms is that particles behave as if they 
were ideal fermions in the box of size L — Naio- In- 
deed, this interpretation explicitly applies to the hard- 
rod gas, where the excluded volume correction is nega- 
tive, as aiD > 0. For a negative scattering length the 
"excluded volume" correction changes its sign and be- 
comes positive L L + N\aiD\- In Fig.[T]we present the 
characteristic behavior of the one-body scattering solu- 
tion at low energy for three short-ranged potentials. The 
Tonks-Girardeau potential corresponds to zero-range in- 
finitely strong repulsion. This places the node of the wave 
function at the origin and according to definition [U the 
value of the s-wave scattering length is zero am = 0. 
For a hard-rod interaction potential the position of the 
first node is positive and thus am > 0. The slope of 
ipif) is determined by the scattering momentum (refer 
to Eq. ^ and the discussion above it), so for a similar 
scattering energy the only relevant difference in the wave 
function corresponding to different interactions is just a 
shift in abscissas. Thus, two Tonks-Girardeau particles 
separated by a distance r and two hard-rod particles sep- 
arated by a distance r — am "feel" each other in the same 
way. The only differences appear at very small distances 
of the order r w aio- In a similar way the scattering so- 
lution for a Lieb-Liniger Hamiltonian can be "adjusted" 
to match the Tonks-Girardeau solution by the change 
r — > r — aiD (mind that aio < in this case). This 
makes natural that the "excluded volume" correction is 
encountered for different potentials and it can change its 
sign. This was first noted in ReL ^4Q] . 



IV. TWO-DIMENSIONAL SYSTEMS 
A. Overview of the equation of state 

A series of our previous works have been devoted to the 
study of the equation of state of dilute two-dimensional 
Bose gases [24 |4ll |42|, |43|, y , |45l El . A number of inter- 
action potentials (dipolar, Yukawa, etc.) were considered 
in a wide range of densities. There it was demonstrated 
that the dipoles crystalize at large densities. The den- 
sity of the quantum phase transition turned out to be 
extremely large nag « 2900 [Uli^,!!!!. This shows that 
dipolar interaction potential is rather "soft" compared to 
the hard-core potentials, which are expected to crystal- 
lize at values of the gas parameter smaller than unity (for 
example, na^ = 0.33(2) in the case of hard-disks '49!|). 
The equation of state of a dilute gas was obtained both 
in the universal and non-universal regimes from Monte 
Carlo calculations for dipoles [4l|, |44[, hard- and soft- 
disks [soj . A peculiarity of the two-dimensional systems 
is that the Gross-Pitaevskii equation has a limited appli- 
cability even in very dilute systems [11] due to the loga- 
rithmic dependence on the gas parameter rather than on 
powers of it, like in three- and one- dimensional systems. 
This means that it is extremely difficult to study numer- 
ically the universal equation of state. Densities as low as 
na^ « 10"^°" had to be reached in calculations [23| in 
order to check numerically the low-density expansion of 
the universal equation of state. It turns out that in order 
to describe correctly the beyond mean-field effects, sev- 
eral terms have to be summed even at densities as low 
as nttQ ~ 10"^" since the series comes out in terms of 
the slow converging logarithm function Inncg, which at 
such densities is of the order of the next (constant) term. 
Historically it turned out to be very difficult to obtain 
the correct expression for this term, see Ref. |2J] for a 
summary of different results. Only recently the correct 
expression was obtained [2l], [S^, [23, |5l[ . 

Once the structure of the universal terms is well es- 
tablished, we are ready to test the concept of an energy- 
dependent s-wave scattering length. 

B. Universal terms 

In this section we study the equation of state of 2D 
Bose gases in the universal regime, i.e. where the in- 
teraction potential can be described by one parameter, 
namely the s-wave scattering length, and all properties of 
the gas are fully defined by the gas parameter na^. With- 
out going through a rigorous derivation of the equation of 
state, which would be an extremely tedious calculation, 
we provide some simple ideas that give insight into the 
relevant physics involved in the equation of state. 

In a weakly interacting system of any dimensional- 
ity the leading contribution to the energy comes from 
the mean-field theory. Assuming that the density is 
low enough, it does not matter what the exact shape 
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of the short-range interaction potential is, and a simple 
(5-pseudopotental can be used. In this way the real in- 
teraction potential can be replaced by a zero-range one, 
such that it imposes a correct zero boundary condition 
to the scattering state 



V{r) — g5{r) x [regularization) 



(9) 



The regularization operator is needed to make the 5- 
function description compatible with a generic 1/r (or 
1/v^) divergence in a 3D (or 2D) geometry, although 
this is not important for our considerations. The sub- 
stitution of ([9]) into the expression of the interaction 
energy written in first quantization simplifies the dou- 
ble integration E ^1/2^ dri J dr2^t(i.^)^t(r2)T/(|ri - 
r2|)^(r2)*(r2) = g/2 Jdr *t(i.)^t(i.)^(i.)^(r). Treat- 
ing the field operator ^(r) as a classical field and sub- 
stituting it with the particle density ^/n one obtains the 
mean-field expression for the energy 
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It is easy to see from Eq. © that the coupling constant 
has dimensionality of [E x L^] and it has to be expressed 
in terms of the parameters of the scattering problem, 
which are h, m and a. In the three-dimensional case the 
considerations of units leads to a combination propor- 
tional to the s-wave scattering length g^D oc h^asD/nT-- 
Indeed, the exact expression is g^u — Airh^a^jj/m. In 
a one-dimensional system the correct units are obtained 
in a combination which is inversely proportional to the 
s-wave scattering length gi^ cx h? / (mauj). This agrees 
with the exact result gm = —2h?/{maiD). The two- 
dimensional case is special in the sense that combina- 
tions having the proper units can be obtained without 
involving the s-wave scattering length g2D oc /m. The 
dependence on a can come only in a combination with the 
scattering momentum k, which in a homogeneous system 
is related to the density. The exact result [52, 53] indeed 
has the anticipated structure g2D = Airh^ / {m\\iina\jj\). 
As explained, the corresponding mean-field term can be 
split into a part independent of a 
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and a part that depends on a 
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The momentum dependence of the 2D coupling con- 
stant makes perturbative theory very involved. More- 
over, the dependence of the coupling constant on the 
s-wave scattering length is very weak as a enters in 
a logarithm. Already at the mean-field level one sees 
that the usual relation between the chemical potential 
^MF _ g^j^^ ^Yie energy per particle E^^^ /N — 1/2 gn 
written as E/N = 1/2 ^ is valid only to the leading 



term. Indeed, direct integration of the chemical poten- 



tial E{N) = fi{N')dN' leads to a different result 
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where T{a,x) = t°'~^e~*dt is the incomplete gamma 
function. The energy expansion in the dilute regime 
na'^ — > can be obtained by doing the integration by 
parts or from the large argument expansion of r(0,a;). 
One finds that this generates a number of terms includ- 
ing a contribution to the BMF constant 
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The difference contain terms that are logarithmically 
small which, anyway, exceeds accuracy of the mean-field 
theory, so that at the MF level there is no legitimate rea- 
son to prefer one expression over the other. On the other 
hand terms of this order are important when one studies 
beyond MF terms. 

The most important beyond mean-field terms were ob- 
tained by V. N. PopovjH in 1972 (see also his book 
He obtained a recursive expression relating the chemical 
potential fi and the density n for a given value of the 
inverse temperature f3 
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where e'^{k) — [h^k"^ /2mY + h^k'^^i/m is the Bogoliubov 
spectrum and Eq is of the order of fi^/mrQ, with tq is 
the range of the interaction potential. We write the last 
relation introducing an unknown coefficient of propor- 
tionality Ci such that Eq = Cih^ /ma? . 

At zero temperature quasiparticle excitations are ab- 
sent and the expression simplifies. By solving Eq. (jlSp 
iteratively one obtains the following expression for the 
chemical potential 
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Lozovik and Yudson used in 1978 diagrammatic tech- 
niques to find a recursive relation which relates the chem- 
ical potential and the density fssj 
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Solving recursively Eq. (fT7)) [see also [28] and Appendix 
B in [54I] one generates the first BMF term In | In na^j 
and the second BMF term proportional to Iutt. 

It should be noted that while in a three-dimensional 
system Bogoliubov approximation works in a dilute 
regime, in the two-dimensional case even in a very rarified 
gas the Bogoliubov theory becomes inapplicable [s^. 

Furthermore, the condensate fraction has a logarithmic 
dependence on the gas parameter hq/u = 1 — 1/| In na^j 
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[5 3 . ISq . The difference between the total density n and 
the condensate density hq in Eq. pU)) leads to an addi- 
tional contribution in the second BMF (constant) term. 

Another contribution to BMF terms should contain 
Euler's 7 constant. The calculation of the chemical po- 
tential summing the ladder diagrams relates fi to the 
scattering amplitude 1 1511 . which itself contains the 7 con- 
stant in the expansion [57[. 

Summarizing, one expects to find the following types 
of BMF corrections: 

• A first BMF term of the form In | In na^ \ 

• A contribution to the second BMF (constant) term 
proportional to In tt 

• A contribution to the second BMF (constant) term 
proportional to 7 

• A contribution to the second BMF (constant) term 
proportional to a constant of the order of 1 

Historically it took a long time to obtain correctly 
BMF expansions at low densities (for a literature review 
refer to [23l)- The double logarithm term can be obtained 
from the iterative relation (jl7p and it is present in the ma- 
jority of theories. Unfortunately, this term alone is not 
sufhcient to describe the universal regime and the cal- 
culation of the all the contributions to the second BMF 
term was a challenging task. We note that the corre- 



sponding problem in 3D was solved in the 1950s [lj| and 
the ID problem in 1960s 

The universal equation of state for the chemical poten- 
tial should read then 
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Notice that this expression is compatible both with 
Eq. (fT6|) and the result of iterating Eq. (fT7|) for ^. The 
second BMF term was recently obtained analytically 
[21I (2^ I2I as Cf = - luTT - 27 - 1 = -3.30... and its 
value was confirmed numerically in Ref . [ISl • The subse- 
quent constant was derived a short time ago in Ref. (5ll | 
with its value given by — —0.751. In the following we 
will use a value obtained from a fit to Monte Carlo data 
C2^ = -0.3(1) [24]. 

The expansion of the energy per particle takes then a 
form similar to (fT8|) 



E 
N 



2Trh^n/m 



I In? 



In I In ) 



In I Inria2| + C| 



with the coefficients related as Cf = Cf + 1/2 = -2.80... 
and Ci = + 1/4: (equals to -0.05(10) from the nu- 
merical fit). 
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Figure 2: (Color online) Finite-energy s-wave scattering 
length as a function of momentum k of the incident particle for 
different interaction potentials. All quantities are measured 
in units of ao. Solid line, hard disks a{k) = ao; dashed lines, 
soft disks: thick line, numerical solution as the node of (|C5|) . 
thin line, analytical expansion a(fc)/ao = 1 — 5.53454A;^ao as 
comes from (|C9|I using the parameters of soft disks taken from 
Ref. 50]; dash-dotted lines, dipoles: thick line, numerical so- 
lution, thin line, fit PU)) . 



C. Non-universal terms 

The specific details of an interaction potential become 
important when the density is large, so that the equa- 
tions of state is no longer universal. In the regime of high 
densities one parameter, namely, the zero-energy s-wave 
scattering length a is no longer sufficient to describe the 
system properties. In Section |TT] we have formulated our 
proposal using an energy dependent s-wave scattering 
length. This allows us to generate non-universal terms 
in an energy expansion and also to understand analyti- 
cally at which densities deviations from the universal law 
appear. This can be applied at densities for which the 
universal equation of state is known. As the reference 
equation of state we take (fTO|) . In this section we test 
our proposal for three different potentials, such as hard 
disks, soft disks and dipoles. 

In the case of hard disks, the interaction has only one 
length scale, namely, the size of the disk. As a result the 
energy dependence is trivial an oik) — oq. The absence 
of additional length scales means that out of all interac- 
tion potentials the hard core potential has an equation 
of state that is the most similar to the universal one. 

In the case of soft disks, corrections due to the fi- 
nite scattering energy are important at typical densi- 
ties na^ > 10~^ [50] ■ The first correction due to the 
finite value of the scattering energy is quadratic in mo- 
mentum, as shown in Appendix [C] The explicit expres- 
sion for a{k) is given by formula (jC9p and it reduces to 
asoik) = asD{0){l - 5.53454^^ + ...) for the choice of 
soft disk parameters as in Ref. [50| . 
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The dipolar interaction potential decays slowly and 
deviations from the universal equation of state appear 



much earlier. In a very dilute system, 



< 10- 



the 



In order to find non-universal corrections to the en- 
ergy, according to the proposed scheme, we substitute 
the gas parameter nog in the universal expansion (fT9|) 



dipole-dipole scattering length is well approximated by with na (k). Within the level of accuracy of interest, 

it is sufficient to use the mean-field expression for the 
scattering momentum cx 2mE /Nh^ — AiTn/\ Inna^j. 



its value at zero scattering momentum, add(O) = e^'^rd — 
3.17222. ..rd, where rd is a characteristic lengthscale for 
dipole-dipole interaction potential [45|. We solve the s- 
wave scattering problem numerically and find that the 
following fit describes well the numerical data for the 
value of the s-wave scattering length at low energies 

add(krd) _ g-cxp{0. 44108240. 31414 In fcrd-0. 0275752 In^ '='"d}(20) 

ad<i(0) 



In the case of soft disks this leads to the substitution 
Inna^ InrtOg -|- 21n(l — a/c^flg). The logarithm can be 
expanded as ln(l — £) = £ + leading to non-universal 
corrections of the order of 2ak'^a'^ oc Snana? /\hina?'\. 
The resulting equation of state for soft disks than reads 



E 
N 



2-Kh^n/m 



|lnna2| +ln|lnna2| -InTr- 27- 1/2+ [ln|lnna2| - Iutt - 27 + 2.0(1) + 1/4 4-87 

I 



(21) 
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Figure 3: (Color online) Energy per particle, analysis of non- 
universal beyond MF corrections. Main figure. Beyond MF 
terms in the energy as a function of the double logarithm of 
the gas parameter. Symbols, DMC results: up triangles, hard 
disks, squares, soft disks, down triangles dipoles. Correspond- 
ing lines: hard disks, Eq. p9|) . soft disks, Eq. (|21|) . dipoles 
Eq. (|19p with a{k) as in Eq. (|2U|I . Inset, energy per particle 
E/N in units of "universal" equation of state, Eq. (|19p . as a 
function of the gas parameter na^. 



is universal shrinks in the case of soft disks and dimin- 
ishes further for dipoles (compare to the dependence of 
the corresponding a(fc). Fig. [2]). We find that the analyt- 
ical description we obtain for the non-universal behavior 
works rather well. In particular, the density, at which de- 
viations from the universal law start to be visible, is pre- 
dicted correctly by our approach. The analytical formula 
(I2ip provides not only a good qualitative description, but 
even the quantitative agreement is good. From Fig. [3] it 
might seem that the description is better for soft disks 
compared to dipoles, but in reality the description for 
both potentials is expected to have a similar level of ac- 
curacy. In order to check that we solved self-consistently 
the equations E = E{n,u,a);a = a{E);u = u{n,a), 
(where u is a dimensionless in-medium scattering am- 
plitude [iJl , see also Appendix |B]) thus obtaining a dif- 
ferent expression, which have the same significative per- 
turbation terms, but differ in higher order terms which 
are outside of the accuracy of our approach. The self- 
consistent solution improves coincidence for dipoles, but 
also changes the predictions for the soft disks introducing 
deviations similar to the ones of dipoles in Fig. [3] 



V. DISCUSSION 



The obtained analytical expressions for the equation 
of state are confronted with the results of Monte Carlo 
simulations. Figure [3] shows the beyond MF energy as a 
function of the double logarithm of the density for differ- 
ent interaction potentials (compare to Fig. 1 in Ref. [2^). 
As anticipated, the beyond MF terms have the most sim- 
ple dependence for hard disk potential (since the s-wave 
scattering length dependence is flat, see Fig. [J) and it 
is the best one described by the "universal" equation of 
state The region where the description of the energy 



Recent progress in techniques of cooling and confine- 
ment permits to realize extremely dilute gases in the 
regime of quantum degeneracy, thus providing a very ad- 
vanced tool for studying properties of weakly interacting 
gases. The s-wave scattering length can be controlled by 
use of Feshbach resonance and can be set to, essentially, 
any desired value by choosing an appropriate magnetic 
field. Many features of the equation of state can be in- 
ferred from measuring energetic properties, such as re- 
lease energies in time of flight experiments. Also the size 
of the cloud and the density proflle are related to the 
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equation of state. The most precise technique for the mo- 
ment is the accurate measurement of the frequencies of 
collective oscillations. This method was successfully used 
to study beyond MF terms in the equation of state of two- 
component Fermi gases in the BCS-BEC crossover [s^. 

In previous sections we have investigated the proper- 
ties of low-dimensional weakly interacting Bose gases as a 
function of the one- and two-dimensional gas parameters 
n-iDaiD and ^21)02^5, respectively. The low-dimensional 
system can be realized in experiments by strongly squeez- 
ing the gas in one or two directions. Assuming that the 
trapping is harmonic with frequency uj the condition of 
being in a low-dimensional regime is that the oscillator 
levels should not be excited neither by the energy per 
particle nor by the temperature E/N, ksT <C fko. 

In a one-dimensional system a relation between the 
three-dimensional a^D and the one-dimensional aiD s- 
wave scattering lengths was found in Ref. [s^ assuming 
harmonic radial confinement with oscillator length aho- 
The relation has a resonant behavior when a^D is of the 
same order as aho due to the contribution of virtual ex- 
citations of the levels of transverse confinement. The 
one-dimensional coupling constant guj = —2h^/{raaiD) 
is expressed as [sof 




Figure 4: (Color online) Comparison of the lowest breathing 
mode frequency, as a function of the coupling strength 
N^^'^ / for different interaction potentials. Solid line, 
hard disks, Eq. (|19|l : dashed line, soft disks, Eq. (121(1 : dash- 
dotted line, dipoles, Eq. (fT9)) with a{k) as in Eq. pOjl. 



1 



giD 



maf^^ 1 - 1.0326 asn/a/io 



(22) 



In particular, at the top of an Olshanii resonance, the 
one-dimensional s-wave scattering vanishes and am — 
0, making gio — > 00. This corresponds to the Tonks- 
Girardeau limit. Close to the resonance is small and 
expansions like (|5l7p are applicable. 

In a similar way to Eq. (j22p . the coupling constant in 
a q uasi-two-dimensional system has a resonant structure 

m 



47rfi2 



92D 



I \n(2Tr\ ij\mal^/h?)\ + V^naho/asD 



(23) 



and describes a competition between "purely two- 
dimensional" logarithmic term and a mean-field Gross- 
Pitaevksii term = 2V27rfi,^/m, which can be ob- 

tained from the Gross-Pitaevskii energy functional as- 
suming a Gaussian profile in the tight direction of the 
confinement. The results from Section IIVI apply to a 
purely two-dimensional system, when the logarithmic 
term in Eq. is dominant, i.e. when a^o ^ o-3D- We 
have explained that the mean-field (here, in a "purely" 
two-dimensional sense) regime is achieved when the dou- 
ble logarithm of the two-dimensional parameter is large 
In I In na^ \ 3> 1 which leads to extremely rarefied densi- 



ties, such as 



< 10 



-862 



Fortunately, the nature pro- 



vides a way to get such small effective densities. Indeed, 
Eq. (j23p can be rewritten introducing the second term 
of the denominator under the logarithm. The resulting 
expression can be interpreted in the sense of a "purely 
two-dimensional" system, but with a rescaled effective 



density n* cx exp{—\/2Traho/cL3D}n. Here the effect of 
large aho/asD ratios is exponentially amplified. 

Expressions for the (quasi)low-dimensional coupling 
constants (|22|23|) were obtained from the analytic solu- 
tion of the two-body scattering problems in the presence 
of a tight harmonic confinement. The existence of one- 
dimensional resonance in a many-body system was later 
confirmed in numerical simulations |6l| . A similar two- 
dimensional study is more involved as the expression of 
the coupling constant depends on the chemical potential 
and we are not aware of such studies. 

The energetic properties of trapped gases can be ac- 
cessed by observing the frequencies of collective oscilla- 
tions. By displacing the center of the trap it is possible to 
generate oscillations which depend only on the frequency 
of the trapping potential. Instead, a sudden change in 
the frequency of the trap causes "breathing" oscillations, 
for which the frequency depends on the compressibility 
of the gas, which, in turn, is related to the equation of 
state. In Fig. |4] we show predictions for different inter- 
action potentials in two-dimensional systems. Out of all 
considered model interactions, the hard core potential 
shows the strongest dependence. The soft disk and dipo- 
lar potentials have softer dependencies. 



VI. CONCLUSIONS 

To conclude, we have studied the energetic proper- 
ties of dilute low-dimensional interacting Bose gases at 
zero temperature. In the regime of ultralow densities 
the equation of state is described by only one dimension- 
less parameter, the gas parameter na^ . The universal 
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equation of state in three- and one- dimensional systems 
dates back to 1960s |1, i, [11. The beyond mean-field 
terms of the two-dimensional equation of state were ob- 
tained recently in Refs. [ll], [S^, [5^ and their correctness 
was verified in Monte Carlo calculations [24l|. When the 
density is increased, the details of the interaction poten- 
tial become important and deviations from the universal 
behavior are observed. 

We propose to use an energy-dependent s-wave scatter- 
ing length to describe the non-universal behavior. This 
method permits to generate non-universal terms in the 
equation of state. The advantage of the proposed ap- 
proach is that it sufficient to known the dependence on 
the s-wave scattering length for this method to be appli- 
cable. This permits us to use it in low dimensional sys- 
tems, where the equation of state in term of non-universal 
are not well established. We test this approach on one- 
dimensional systems, where a direct comparison to ex- 
actly solvable models is done, and on two-dimensional 
systems, where numerical results for different interaction 
potentials (hard disks, soft disks, dipoles are used. We 
find that the typical density at which the non-universal 
terms become important is correctly estimated. For one- 
dimensional systems, in the cases when the energy ex- 
pansion can be obtained exactly, we show that the struc- 
ture of the potential-dependent terms is predicted cor- 
rectly. Finally we point out that non-universal terms 
can be studied experimentally by observing frequencies 
of collective oscillations. 
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equations 

e(7) 



^ 1 



7 = 




p(fc) 



2\p{x) dx. 



2tt J A2 + (fc - 2n ' 

-1 



(Al) 
(A2) 
(A3) 



where 7 = —2/nijjaiD- 

It is possible to obtain explicit expressions for the en- 
ergy in terms of the gas parameter in the limits of small 
and large gas parameter as a series expansion. We solve 
the system of Eqs. (|AmA3[) iteratively in the regime 
\na\ <C 1. This is done by starting from — l/(27r) 
and substituting it into the r.h.s. of Eq. (jASp to get p^^^^ 
The obtained expression is used for the next iteration 
and so on. 

We provide an explicit expression for the ground-state 
energy close to the Tonks- Girardeau regime: 



E 



6m 



.i=0 



327r^ ^, _4, 
+ 1^ + 0(7 1 



(A4) 



The first two terms: E/N ^ 



4/7) were 



obtained in the original work of Licb and Liniger [J] (see 
also [Hil]). 

We also rewrite expression (|A4p in terms of the gas 
parameter h = nnyaiD as 



E 

N 



6m 



47r"^a 



2 ~ 3 



Li=o 



15 



0(^4 



(A5) 



The "excluded volume" contribution is intentionally sep- 
arated from the non- universal part. 



Appendix B: EQUATION OF STATE OF A 
TWO-DIMENSIONAL BOSE GAS FROM 
CHERNY AND SHANENKO THEORY. 



Appendix A: EXACT EQUATION OF STATE OF 
WEAKLY-INTERACTING ONE-DIMENSIONAL 
BOSE GAS WITH 5-PSEUDOPOTENTIAL 
INTERACTIONS. 



We note that the derivation of the equation of state of 
a weakly interacting Bose gas, proposed by Gherny and 
Shanenko [2l|, can be used to obtain an explicit expres- 
sion of the energy as a function of the density n: 



E _ 2TTh^n{u) 
N m 



u 
~2 



2v?e 



2/u 



Ei - 



n{u)a^ 



exp{— 27 — l/u} 



(Bl) 
(B2) 



The ground-state energy of a Lieb-Liniger gas can be 
found exactly using the Bethe ansatz approach. The 
energy as a function of the gas parameter is obtained 
implicitly by solving the following system of integral 



here u (dimensionless in-medium scattering amplitude) 
defines the parametrical dependence of the energy on the 



density, Ei(a;) 



function, and 7 is Euler's constant. 



e */t dt is exponential integral 
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Appendix C: FINITE-ENERGY SCATTERING 
PROBLEM FOR A SOFT DISK POTENTIAL IN 

2D. 

In order to find the two-dimensional s-wave scattering 
length of the soft-disk potential in 2D we have to solve 
the two-body scattering problem. The positive energy 
Scrodinger equation for two particles of equal mass reads 

--A/(r) + F(r)/(r) = ^/(r), (CI) 
m m 

where k is the relative momentum. We consider the soft 
disk interaction potential and look for a spherically sym- 
metric solution. The interaction potential is defined by 
the range of the potential Rq and the height of the soft 
disk by 



where k = — k"^. The s-wave scattering length is the 
node of the function (|C5p closest to the origin. We will 
consider the case of low densities, so that the incident 
particles is slow, k <^ x. Then /(r) can be expanded in 
powers of k and one has 



+ ^[(l-rVi?^)/o(-i?o) 

2 + K^iR^o - r2) - x2(r2 + R^) ln(i?o/r) 



>t:Rn 



Oik') 



hixRo) 
(C6) 



V{r) 



h^x^/m, \r\ < Ro 



0, 



\r\ > i?o 



(C2) 



The zero energy s-wave scattering length is found by set- 
ting fc = in (jC6p and leads to 



In the inner region, r < i?o, we use a solution that is 
regular at the origin, 



fir) = /o(rV^2_fc2)^ < 



(C3) 



where /o(a;) is the modified Bessel function of the first 
kind. The normalization constant is not important for 
the present considerations. In the outer region the solu- 
tion is simply a two-dimensional plane wave 

fir) - CiJoikr) + CWoikr), \r\ > Ro, (C4) 

where Joix) and lo(x) are Bessel functions of the first 
and second kind, respectively. The coefficients Ci and 
C2 are obtained from the continuity condition for /(r) 
and fir) at the edge of the soft disk r — Rq. This gives 
the following solution in the outer region, r > Rq. 



fir) 



ttRq 



{KhinRo) [JoikRo)Yoikr) ~ Mkr)YoikRo)] 



khinRo) [JiikRo)Yoikr) - Joikr)YiikRo)]} ,(C5) 



hq = exp 



Jo(^i?o) 



Ro 



(C7) 



xRohixRq) 

Furthermore, one can set r — > ao in the second line of 
Eq. (|C6p and find a correction to the position of the node 



aik) = ao- ak'^al + Oik'^), (C8) 

Ro loi^Ro) + hixRo) _ 1 

4 



hi^Ro) 



Ro 



1 



(C9) 



In order to test the accuracy of Eq. (|C9p we also find nu- 
merically the nodes of the finite-energy scattering func- 
tion (|C5[) . A comparison of the exact result for the s-wave 
scattering length a(A;) to the expansion (|C9[) is presented 
in Fig. [21 We see that at the densities of interest, the 
obtained expansion works very well. 
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